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ABSTRACT 

Aims. We report the first detailed X-ray study of the supernova remnant (SNR) G304.6+0. 1, achieved with the XMM-Newton mission. 
Methods. The powerful imaging capability of XMM-Newton was used to study the X-ray characteristics of the remnant at different 
energy ranges. The X-ray morphology and spectral properties were analyzed. In addittion, radio and mid-infrared data obtained with 
the Molonglo Observatory Synthesis Telescope and the Spitzer Space Telescope were used to study the association with the detected 
X-ray emission and to understand the structure of the SNR at differents wavelengths. 

Results. The SNR shows an extended and arc-like internal structure in the X-ray band with out a compact point-like source inside 
the remnant. We find a high column density of Njj in the range 2.5-3.5xl0 22 cirT 2 , which supports a relatively distant location (d > 
9.7 kpc). The X-ray spectrum exhibits at least three emission lines, indicating that the X-ray emission has a thin thermal plasma 
origin, although a non-thermal contribution cannot be discarded. The spectra of three different regions (north, center and south) 
are well represented by a combination of a non-equilibrium ionization (PSHOCK) and a power-law (PL) model. The mid-infrared 
observations show a bright filamentary structure along the north-south direction coincident with the NW radio shell. This suggests 
that Kes 17 is propagating in a non-uniform environment with high density and that the shock front is interacting with several adjacent 
massive molecular clouds. The good correspondence of radio and mid-infrared emissions suggests that the filamentary features are 
caused by shock compression. The X-ray characteristics and well-known radio parameters indicate that G304.6+0. 1 is a middle-aged 
SNR (2.8-6.4)xl0 4 yr old and a new member of the recently proposed group of mixed-morphology SNRs. 

Key words. ISM: individual objects: G304.6+0.1 - ISM: supernova remnants - X-rays: ISM - Radiation mechanisms: thermal 



• tH . 1. Introduction Observatory Synthesis Telescope (MOST) Catalog (Whiteoak & 
/\ Green 1996) as a small-diameter SNR with irregular shell mor- 
»H ; The supernova remnant G304.6+0.1 (also known as Kes 17) was phology. In addittion, Frail et al. (1996) have detected 1720 MHz 
. TO , first detected at 408 MHz and 5 GHz by Shaver & Goss (1970), OH maser emission around the object using the Green Bank tele- 
using the Parkes and Molonglo radio telescopes. The authors ob- SCO p e 
tained spectral information and some physical characteristics of 
the SNR. A spectral index a 0.51 (S oc V"), an angular size 

of 8.4 arcmin and a possible distance of 12 kpc were estimated More recently, Reach et al. (2006), using mid-infrared data 
for the source. Afterwards, Milne & Dickel (1975) determined from me Spitzer Space Telescope (Werner et al. 2004) at 3.6, 
a radio polarization of 1.3%. Using the Parkes hydrogen line 4,5, 5 >8> and 8/Um of G304.6+0.1, revealed bright thin isolated 
interferometer, Caswell et al. (1975a) constrained the distance n i ame nts in each of the IRAC channels. Based on colors and 
to a minimum of 9.7 kpc. Later, Milne et al. (1985) mapped the detailed morphological agreement of the images in the four 
the object with the Fleurs array, and suggested that its peculiar channels, the authors suggested that most of the mid-infrared 
structure resulted from the shock front collision with a dense in- emission from the shell could be produced in the shocked molec- 
sterstellar cloud. The source was then reported in the Molonglo ular mate rial. Recently, Hewitt et al. (2009) computed the total 
hydrogen column density along the line of sight to the SNR, us- 
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Fig. 1. Left panel: XMM-Newton color image of the X-ray emission in the 0.3-10 keV energy range (in blue) of G304. 6+0.1 
with the radio contours (in white) at 843 MHz overlaid. The extraction regions used for the spectral analysis are also indicated as 
yellow circles. Central panel: IRAC three-channels composed image of G304.6+0. 1 with the radio contours (in white) at 843 MHz 
overlaid. Right panel: Composite color image of the radio, mid-infrared, and X-ray emissions of G304.6+0.1. All images are in 
equatorial coordinates J2000. 



interior X-ray emission. However, at present no X-ray study of 
the source has been published. 

We report an X-ray study of SNR G304.6+0.1, using the 
greatly enhanced sensitivity of the XMM-Newton mission. An 
instrument with these characteristics offers a unique opportunity 
to detect X-ray emission from distant and obscured SNRs. Radio 
and infrared observations were also used for a multiwavelength 
study of the source. This work is part of a program aimed to iden- 
tify and study X-ray emission from SNRs. The structure of the 
paper is as follows: in Sect. [2] we describe the XMM-Newton ob- 
servations and data reduction. X-ray analysis and radio/infrared 
results are presented in Sects 3 and 4, respectively. In Sect. |5j 
we discuss the implications of our results. Finally, we summa- 
rize the main conclusions in Sect. 6. 



3. X-ray analysis of G304.6+0.1 

3.1. Image 

In Fig. 1 (left panel) we show the X-ray image of G304.6+0. 1 in 
the 0.5-0.8 keV energy range. The total unabsorbed X-ray flux 
for this band (after the fit to the X-ray spectrum, see Table 1 for 
details) is F X (o.5-8.0 keV)= 7.1 1(+0.36) xl0~ n erg cirT 2 s -1 . 

The image reveals diffuse and extended X-ray emission with 
an irregular filled-center and possible arc-like features. No com- 
pact point-like X-ray source was found within the diffuse emis- 
sion. Assuming a central source with a typical compact central 
object (CCO) spectrum (r ~ 5), we obtain a flux upper limit of 
6.5xl0~ 15 erg cm" 2 s -1 in the 0.5-8 keV band. Clearly the X- 
ray emission seems to correlate well with the internal part of the 
NW radio shell. An irregular region of X-ray emission is also 
observed at (a, (5) J2 ooo.o = (13 h 05 m 35 s .0, -62°44'00"). 



2. Observations and data reduction 

The field of SNR G304.6+0.1 was observed on 2005 August 
25 by the XMM-Newton X-ray satellite during revolution 
1046 (Obs-Id. 0303100201), with a total integration time of 
26.8 ksec. The observation was centered on (a, (5)j2ooo.o = 
(13 h 05 m 46 s .0,-62°43'40"), and acquired by the EPIC-MOS 
(Turner et al., 2001) and EPIC-PN (Struder et al. 2001). The 
data were taken with a medium filter in full-frame (FF) imaging 
mode. The data were analyzed with the XMM Science Analysis 
System (SAS) version 9.0.0 and the latest calibrations. 

To exclude high background flares, which could eventually 
affect the observations, we extracted light curves of photons 
above 10 keV for the entire field-of-view of the EPIC, and dis- 
carded time intervals in which background flares occurred for 
subsequent analysis. We finally used a total of 21.181 ksec, 
23.962 ksec, and 23.957 ksec for the PN, MOS-1 and MOS- 
2, respectively. To create images, spectra, and light curves, we 
selected events with FLAG=0, and PATTERNS = 12 and 4 for 
MOS and PN cameras, respectively. Hereafter we use clean 
event-files in the [0.3-10] keV energy band. 



3.2. Spectral analysis 

The X-ray spectrum of G304.6+0. 1 was extracted from the EPIC 
cameras. We used the SAS task evselect with appropiate pa- 
rameters for the PN and MOS 1/2 cameras. The extraction re- 
gions used were circles with radii 1.4, 1.8, and 0.75 arcmin (see 
Fig. 1, left panel). Background spectra were extracted from a re- 
gion where no X-ray emission was detected. 

Figure 2 shows the background-substracted EPIC-PN 
spectra obtained from the XMM-Newton observation for three 
differents regions (i.e., north, center, and south). We grouped 
the extracted EPIC PN and MOS spectra with a minimum of 25 
and 16 counts per spectral bin. Ancillary response files (ARFs) 
and redistribution matrix files (RMFs) were calculated. At high 
X-ray energies (above 7 keV), the full spectrum shows features 
with low statistical significance (1 to 1.5 cr), which are probably 
related to fluorescence lines in the background spectrum of 
XMM-Newton (e.g. De Luca & Molendi, 2004). The spectral 
analysis was performed with the XSPEC package (Arnaud, 
1996). The X-ray spectrum of the remnant exhibits emission 
lines at energies of 1.54 keV (Fe XXIII+Fe XXIV+Mg XIV), 
1.85 keV (Si XIII), and 2.45 keV (S XV). However, the quality 
of the data prevents a further refined analysis of the abundances 
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Fig. 2. Combined XMM-Newton EPIC PN and MOS X-ray 
spectra for the regions north, center, and south of G304. 6+0.1 
and residuals (botton panel in each spectrum). The solid line in- 
dicates the best-fit pshock+PL model (see Table 1). The iden- 
tified lines, mostly blended, correspond to atomic transitions of 
Fe XXIII+Fe XXIV+Mg XIV, Si XIII, and S XV. 



of each single ion. 

The X-ray spectrum was fitted with various models: (i) a 
pshock non-equilibrium ionization model with constant temper- 
ature (Mazzotta et al.1998); (ii) a mekal emission model from 
hot diffuse gas (Mewe et al. 1985), and (Hi) a single power- 
law (PL) model. Each model was modified by a low-energy ab- 
sorption model WABS (Morrison & McCammon 1983). Several 
combinations were tentatively used to fit the spectra in these 
regions. We initially varied individual parameters: absorption 
(Nh), temperature (kT), F index, ionization time scale (t), and 



normalization during the minimization process. As a first step, 
we tried the spectral fit by using single thermal models, but 
they did not yield acceptable fits at hard energies (i.e. E > 4 
keV), leading to large residual statistics, i.e. xl ~ 2.0 or even 
larger. Otherwise, the non-thermal emission model failed to fit 
the spectrum at the soft energies, affecting the Nh determina- 
tion. We decided to treat the abundance as a single parameter 
in the pshock model. We initially used the solar abundances of 
Anders & Grevesse (1989), which were also thawed during the 
minimization procedure. The goodness of the model fit was de- 
rived according to the 2-test statistics^ 

A combination of two thermal models works well, but the 
temperature of the second (hard) component is mismatched, 
therefore it does not provide a satisfactory description of the 
X-ray spectrum. It reaches temperatures higher than 10 keV, 
which cannot be considered a reasonable solution for a restricted 
0.3-8.0 keV X-ray spectrum. The results suggest that the ther- 
mal emission dominates at low and intermediate energies, while 
a small contribution of non-thermal emission (~8%) is impor- 
tant at high energies. Therefore, on the basis of x 2 minimiza- 
tion, a combination of thermal and non-thermal models provides 
the best acceptable fit description of the X-ray spectra of the 
of G344. 6+0.1 in the three regions. The power-law model ad- 
equately fits high-energy photons (up to 5 keV) in which par- 
ticle acceleration in shock-fronts could physically explain this 
X-ray emission. The parameter values of the thermal and non- 
thermal models are physically reasonable (see Table 1 for de- 
tails). Unfortunately, the sensitivity and spatial/spectral resolu- 
tion is not suitable to isolate the thermal and non-themal X-ray 
contributions. 



4. Radio and infrared characteristics of G304.6+0.1 

To investigate the correlation between the X-ray, radio, and mid- 
infrared emissions of G304. 6+0.1, we used radio observations 
performed at 0.843 GHz (with a resolution of 43") by the MOST, 
and mid-infrared data of the Infrared Array Camera (IRAQ 
(Fazio et al. 2004), obtained by the Spitzer Space Telescope at 
3.6, 4.5, 5.8, and 8/rni, respectively. Figure 1 shows three im- 
ages that allow us to compare the correlation among the radio, 
infrared, and X-ray emissions. Figure 1, left panel, shows the X- 
ray emission of the SNR in the 0.3-10 keV energy range, with 
the 843 MHz radio contours (in white) overlaid. Figure 1, cen- 
tral panel, shows the Spitzer/IRAC mid-infrared color image of 
the region where G304.6+0. 1 lies, with the radio contours at 843 
MHz (in white) superimposed. Finally, a composite color image 
of the radio/mid-infrared and X-ray emission of G304. 6+0.1 is 
shown in Fig.l, right panel. 

At radio frequencies, the source displays an asymmetric bi- 
lateral morphology with clear shells in the northwest, and south- 
ern parts of the remnant. This morphology could be the effect 
of a non-uniform ISM, or of a non-uniform ambient magnetic 
field (Orlando et al. 2007). The remnant has an angular size of 
~ 8x8 arcmin, which for a distance of 9.7 kpc (Caswell et al. 
1975b) corresponds to a diameter of ~ 20 pc. The MOST radio 
flux density at 843 MHz is Ss43mhz=18 Jy, which combined with 
that at 408 MHz S 40 8mhz=29.8 Jy (Shaver & Goss 1970) gives a 
spectral index a 0.49 (S oc V"). 

At the mid-infrared part of the spectrum (Fig.l, central 
panel), we note a bright filamentary structure running along 
the north-south direction on the western side of G304. 6+0.1, 
which correlates very well with the NW radio shell. This feature, 



C-statistics was also applied to the spectral fit with similar results. 
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Table 1. X-ray spectral parameters of G344.6+0.1 



Model 


South 


Center 


North 


Parameters 








WABS 








N H [cm- 2 ] 


2.5(±0.2)xl0 22 


3.2(±0.2)xl0 22 


3.5(±0.3)xl0 22 


PSHOCK 








kT [keV] 


0.85(±0.1) 


0.72(±0.1) 


0.68(±0.09) 


Abundance 


U.J(± U.l) 


n o/_i_n i\ 


l.z(±U.O) 


t[s cirT 3 ] 


3.1(±0.7)xl0 12 


2.0(±0.6)xl0 12 


1.2(±0.7)xl0 12 


Norm 


9.0(±0.2)xl0- 4 


1.3(±0.1)xl0- 3 


6.0(±0.2)xl0- 3 


EM (xlO 57 ) 


5.8327562 


12.142146 


9.1539998 


POWER LAW 








r 


2.4(± 0.6) 


1.8(+ 0.7) 


3.1(+ 0.3) 


Norm 


l.l(±0.9)xl0- 4 


2.0(±0.7)xl0~ 4 


6.1(±0.7)xl0- 4 


^/d.o.f. 


0.76/2030 


0.83/3016 


0.8/2131 


Flux [erg cm 2 s '] 


0.5(±0.2)xl0-" 


3.18(±0.07)xl0- u 


3.43(±0.09)xl0- n 



Normalization is defined as 10 14 /4/rD 2 x J n H n e dV, where D is distance in [cm], n H is the hydrogen density [cm 3 ], n e is the electron density 
[cm -3 ], and V is the volume [cm 3 ]. The flux is absorption-corrected in the 0.5-8.0 keV energy range. Values in parentheses are the single parameter 
90% confidence interval. The abundance parameter is given relative to the solar values of Anders & Grevesse (1989). 



first noted by Reach et al.(2006), is centered at (a, (5)j2ooo o = 
(13 h 05 m 3O!0,-62°41'09") and has an angular size of ~ 6 ar- 
cmin. The extended mid-infrared fluxes of this structure at 3.6, 
4.5, 5.8, and 8yum, after the subtraction of intense point-like 
sources in the region, are: F^^ m ~ 1.4 Jy, F4,^ m ~ 1.46 Jy, 
Fs.^m ~ 7.7 Jy, F$x) Mm ~ 21.6 Jy, respectively. Evidently the 
Western radio peak agrees very well with the central part of the 
main mid-infrared arc-like structure. Another remarkable set of 
small and thin like-arc structure is also clearly noticeable just at 
the inner boundary of the NW radio shell. All these features are 
discernible in the four IRAC channels. 

On the contrary, the diffuse X-ray emission observed by the 
XMM-Newton telescope (plotted in blue) is characterized by an 
irregular centrally filled region, which seems to follow the in- 
ternal shell-like radio/mid-infrared morphology (see Fig.l, right 
panel). Another distinct component of the X-ray emission is a 
faint protruding structure that is clearly revealed at the SW side 
of the radio shell. This X-ray feature is coincident with a ra- 
dio feature with a similar morphology, which probably orig- 
inated in a relatively low-density region of the ISM. Finally, 
note that within this region the bright infrared source IRAS 
13024-6227 is embedded (Kleinmann et al. 1986), located at 
(a,5)j2000.o = (13 h 05 m 34 s :0, -62°43'35"). 



5. Discussion 

The main radio and X-ray characteristics of G304.6+0. 1 are a 
centrally filled X-ray morphology dominated by thermal X-ray 
emission (92 % of the total X-ray flux), which is not well cor- 
related with the non-thermal radio shell. This suggests that the 
SNR belongs to the class of mixed-morphology (MM) remnants 
(Rho & Petre 1996). The non-thermal emission represents only 
a small fraction of the X-ray emission, indicating that only a mi- 
nor part of the flux could be associated with highly energetic 
electrons accelerated at the SNR front shock. The detection of 
OH masers at 1720 MHz by Frail et al.(1996) is an indicator of 
shock interaction with a molecular cloud, a result that supports 
the notion that G304.6+0.1 is a member of the proposed group 
of MM SNRs. A list containing SNRs with these characteristics 
can be found in Rho & Petre (1998). 



There are about 20 SNRs with convincing evidence for in- 
teraction with ambient molecular clouds. A substantial fraction 
of them are MM or thermal composite SNRs, with central filled 
thermal X-ray emission surrounded by radio and (in some cases) 
infrared shell structures (see e.g., Yusef-Zadeh et al. 2003; Yang 
& Yang 2007). When a SN begins, the energetic impact of the 
explosion can drive an asymmetric expanding bubble. If the 
surrounding ambient gas is non-uniform, the shock-front could 
severely affect the dynamic structure of the gas and trigger the 
irregular radio and infrared structures. From Fig.l (right panel) 
it is evident that there are many similarities between the ra- 
dio and mid-infrared morphologies observed on the NW part 
of G304.6+0.1. Additionally, the X-ray emission is anticorre- 
lated with both mid-infrared and radio brightness. These results 
suggest that Kes 17 is propagating on a non-uniform environ- 
ment and that the shock front is interacting with adjacent mas- 
sive molecular clouds. The shell-like structures may be caused 
by shock compression. 

At least four possible scenarios have been introduced in past 
years to explain thermal X-ray radiation inside radio shells of 
SNRs: i) cloudlet evaporation in the SNR interior (White & 
Long 1991), ii) thermal conduction smoothing out the temper- 
ature gradient across the SNR and enhancing the central density 
(Cox et al. 1999), iii) a radiatively cooled rim with a hot interior 
(Harms et al. 1997), and iv) possible collisions with molecular 
clouds (Safi-Harb et al. 2005). 

The source G304.6+0.1 seems to propagate into a region 
of greatly varying density, possibly close to several molecular 
clouds. The detection of 1720 MHz OH maser emission (see 
Frail et al. 1996; McDonnell et al. 2008) near G304.6+0.1 is 
an indicator of the interaction with a molecular cloud. Then 
the most suitable scenario for describing the SNR evolution is 
that developed by White & Long (1991). In these remnants, the 
clouds were overrun by the shock and are now cooling and evap- 
orating by saturated conduction in the hot post-shock region 
(Cowie & McKee, 1977). For this reason, the observed X-ray 
morphology may be very different than that of ordinary shell- 
like SNRs. 

It is possible to compute some physical parameters of 
G304.6+0. 1 , using the radio/mid-infrared and X-ray information 
gathered here. From the X-ray image, we can estimate the vol- 
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ume V of the X-ray emitting plasma for the differents regions. 
Assuming that the plasma fills regions like those observed in Fig. 
1 (circles) and assuming a distance of 9.7 kpc, we obtain vol- 
umes V n orth= 9.5X10 57 Cm 3 , Vcenter= 1.5xl0 58 CIT1 3 , and V SO uth = 

l.lxlO 57 cm 3 , respectively. Based on the emission measure 
(EM) determined by the spectral fitting (see Table 1), we can 
estimate the electron density of the plasma, n e , by n e — y/EM/V, 
which results in « e (north)~ 0.99 cm -3 , n e (centre)~ 0.89 cm -3 , 
and n e (south)~ 2.26 cirT 3 , respectively. In this case, the number 
density of the nucleons was simply assumed to be the same as 
that of electrons. The age t was then determined from the ioniza- 
tion timescale, t, by t-T/n e . Therefore, the elapsed time after the 
plasma was heated is in the range of (2.8-6.4)xl0 4 yr. This result 
shows that G304. 6+0.1 is a middle-aged SNR. The total mass of 
the plasma M tota i was estimated by M tot!i i=n e V mu ~ 8-15 M Q , 
where toh is the mass of a hydrogen atom. 

Finally, we compared the X-ray luminosity of G304. 6+0.1 
with other prototypical MM SNRs. The unabsorbed X-ray lu- 
minosity in the 0.5-8.0 keV band of G304. 6+0.1 for a distance 
of 9.7 kpc is Lx ~ 7.5 x 10 35 erg s _1 . This value is similar to 
the SNR G359.1-0.5 (Yusef-Zadeh et al. 1995), but one order of 
magnitude higher than that of W28 (Claussen et al. 1997), which 
was computed at the same energy range. 



6. Conclusions 

We presented the first detailed multiwavelength study of the 
SNR G304.6+0.1, using radio, mid-infrared, and X-ray observa- 
tions. The X-ray morphology reveals diffuse interior X-ray emis- 
sion with an irregular filled-center and possible arc-type struc- 
ture, which seems to follow the internal part of the NW radio 
shell. The XMM-Newton EPIC spectrum displays several emis- 
sion lines, which support the scenario of a thermal nature for the 
X-ray emission. However, a single-component model fails to de- 
scribe the global X-ray spectrum, and at least two components 
are required. The spectrum is well represented by a combination 
of a thermal and non-thermal models. The mid-infrared observa- 
tions show a bright filamentary structure along the north-south 
direction coincident with the NW radio shell. This suggests that 
Kes 17 is propagating on a nonuniform environment with high 
density and the front shock is interacting with adjacent massive 
molecular clouds. The good correspondence of radio and mid- 
infrared emissions implies that these shells are caused by shock 
compression. 

The radio, mid-infrared, and X-ray characteristics of 
G304. 6+0.1 (centrally filled X-ray morphology dominated by 
thermal X-ray emission, not correlated with the non-thermal ra- 
dio shell) indicate that the object is a middle-aged SNR, and that 
it belongs to the class of mixed-morphology remnants. Future 
Chandra observations could help to reveal more details of the X- 
ray structure and to discern the probable presence of a compact 
point-like X-ray source inside the remnant. 
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